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Abstract 

We consider tripartite entangled states for continuous variable systems of EPR type, which 
generalise the famous bipartite CV EPR states (eigenvectors of conjugate choices Xi —X2, Pi + 
P2, of the systems' relative position and total momentum variables). We give the regularised 
forms of such tripartite EPR states in second-quantised formulation, and derive their Wigner 
functions. This is directly compared with the established NOPA-like states from quantum op- 
tics. Whereas the multipartite entangled states of NOPA type have singular Wigner functions 
in the limit of large squeezing, r 00, or tanhr — > 1~ (approaching the EPR states in the 
bipartite case), our regularised tripartite EPR states show singular behaviour not only in the 
approach to the EPR-type region (s ^ 1 in our notation), but also for an additional, auxiliary 
regime of the regulator (s — > \/2)- While the s — > 1 limit pertains to tripartite CV states with 
singular eigenstates of the relative coordinates and remaining squeezed in the total momen- 
tum, the s — > \/2 limit yields singular eigenstates of the total momentum, but squeezed in the 
relative coordinates. Regarded as expectation values of displaced parity measurements, the 
tripartite Wigner functions provide the ingredients for generalised CHSH inequalities. Viola- 
tions of the tripartite CHSH bound {B^ < 2) are established, with B^ = 2.09 in the canonical 
regime {s — > l"*"), as well as B^ = 2.32 in the auxiliary regime {s — > \/2+). 
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1 Introduction 



The nature of quantum entanglement has been pursued almost since the inception of quantum 
mechanics itself. Whereas the early insights of Schrodinger, as well as of Einstein, Podolsky 
and Rosen (EPR) [1] were framed in the Gedankenexperiment mode of discussion with con- 
tinuous degrees of freedom (particle position and momentum eigenstates), the issues were 
taken up quantitatively in Bell's theorem [2| for the case of spin degrees of freedom, via the 
transcription to this context given by Bohm [3]. However, continuous variable (CV) systems 
are the natural framework for most quantum optics and quantum communication work [U [5] , 
and Bell, and the more general Clauser, Horne, Shimony and Holt (CHSH) inequalities [6], 
are important measures of entanglement. A technical difficulty in working with continuous 
variables is that the theoretical ideal EPR type states are singular, whereas experimental in- 
vestigations require regularised states. In the bipartite case these may be provided by so-called 
NOPA squeezed states [3 El El lOl [H] in the large squeezing limit. In an alternative approach, 
Fan and Klauder [T2] constructed somewhat more general classes of EPR states, but without 
providing a regularisation. 

In the extension to multipartite cases, a natural question is the choice of relative variables, 
chosen from amongst the positions and momenta of the constituent particles of the system, 
which will provide an appropriate generalisation of the bipartite EPR states (with or without 
regularisationjl. One possibility is provided by the multipartite Greenberger-Horne-Zeilinger 
NOPA-like squeezed states [131 [H]. These have the virtue of experimental accessibility [8l[9], 
and do show singular behaviour in the large squeezing limit which moreover leads to violations 
of the multipartite CHSH inequalities [l5j . For a derivation of the CHSH inequalities for A^- 
particle systems see for example |16] . However, many other choices of relative variables exist 
- see for example [17] and references therein. 

In this paper we take up a logical generalisation of the original EPR suggestion, in selecting 
simultaneously diagonalisable joint degrees of freedom from amongst the canonical Jacobi 
relative coordinates of the particles. §2 is divided into two subsections. In the first we introduce 
the bipartite Fan and Klauder EPR-like state and propose a possible regularisation. This is 
shown to be identical to the bipartite NOPA state (which approximates the ideal EPR limit) 
with squeezing parameter r — > oo. In fact the correspondence is that our regularisation 
s — > 1"*" coincides with tanhr — > 1~ as r — > oo, with tanhr = In the bipartite case the 

Fan and Klauder states are more general than the original EPR states in that they realise 
explicit nonzero eigenvalues of total or relative positions and momenta; but a study of the 
Wigner functions reveals that such nonzero eigenvalues can be absorbed into shifts of the 
complex displacement parameters, and so in the bipartite case the NOPA states do not lose 
any generality in not allowing for such nonzero eigenvectors. In the next subsection we follow 
the Fan-Klauder approach, in second quantised formalism, to derive the explicit theoretical 
tripartite CV states of EPR type conforming to this structure, and we develop the methods 
to provide a plausible regularisation. 

In §3, we derive explicit Wigner functions for our regularised tripartite states of EPR type 
by interpreting the Wigner functions themselves as expectation values of displaced parity 
measurement operators. As could be expected from their different second-quantised forms, 
the tripartite EPR and NOPA Wigner functions differ significantly (the appendix provides a 
comparison of the second-quantised forms and their respective Wigner functions, including an 
explicit evaluation of the former in the NOPA-like case, and a detailed derivation of the latter 
for our EPR- type states). Specifically, whereas the multipartite NOPA Wigner functions are 
singular in the large squeezing limit, our regularised tripartite states of EPR type admit two 
different singular regimes: not only in the EPR-type regime (s ^ 1 in our notation), where 

■'We refer conventionally to the subsystems as 'particles', but it should be borne in mind that the CV 
systems could equally be independent photon polarisation modes, photon modes or even joint photon and 
phonon degrees of freedom. 
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of course the Wigner function still differs from that of NOPA, but also, rather unexpectedly, 
for an additional, auxiliary regime of the regulator (s V^). In §4 we exploit the fact that 
Wigner functions are immediately applicable as summands in the appropriate tripartite CHSH 
inequalities. We explore the two singular regimes and their Bell operator expectation values 
which control the classical-quantum boundary via the CHSH bound {B3 < 2) and identify 
some instances of violations for each of the cases s — > 1^ and s — > \/2+. §5 includes further 
discussion of our findings, as well as comparison with the recent work [18', '19], which provides 
a general construction of ideal EPR states, and some concluding remarks. 



2 Regularised CV EPR states 
2.1 Bipartite states 

The case considered by EPR in [1] discusses the simultaneous diagonalisation of the two 
commuting variables of difference in position (Xi —X2) and total momentum (Pi + ^2)1 where 
Xj, Pj, j = 1,2 are a standard pair of canonically conjugate variables with [Xj,Pfe] = i6jk- 
Fan and Klauder [12] give an explicit form for the common eigenvectors of the relative position 
and total momentum for two EPR particles in terms of creation operators as follows: 

\rj) =e-|H'+'''*^-''*^^+'^^*^|00), (1) 

where rj = rji + ir/2 is an arbitrary complex number, [a, a^] = 1, [6,6t] = 1 and |00) = |0,0), 
the two-mode vacuum state. Thus 

(Xi - X2) |r?) = V2m\v), (Pi + P2) |r/) = ^2^]^). (2) 
with the coordinate and momentum operators definable as: 

As a genuine representation of ideal generalised EPR states, with appropriate orthonor- 
mality and completeness, \rf) is singular. In this paper we consider the following regularised 
version: 

Ir,), := iV2e-^l^l'+^''"'-^^''*'''+^"'''>0), (4) 

with normalisation \N2\^ = — l)|^/^/s^. 

The bipartite CV state @ is to be compared with a regularised EPR-like state which has 
already appeared in the literature - the so-called NOPA state from quantum optics [7]-|llj: 

\NOPA) = eK'^^^^-«'')|00). (5) 

NOPA states are produced by Nondegenerate Optical Parametric Amplification, and are the 
optical analog to the EPR state in the limit of strong squeezing. The NOPA state has already 
been shown to be a genuinely entangled state that produces violations of the CHSH inequality 

[iniEiiE]. 

Following [50] on reordering SU{i, 1) operators, we can reorder the expression for NOPA 
([5]) into the following form: 

\NOPA) = e''(''^*^-°'')|00) 

_ gratfet^_21n cosh(r)i(ata+fetfe+l)^— rafejQg^ 

= Vl-tanhV e*^"^^™^''^|00). (6) 
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Note from ([6]) that in the number basis, as tanhr — > 1 the \NOPA) state approximates the 
ideal EPR Umit: 

Um \NOPA) « \EPR) « |0,0) + |1,1) + |2,2) + .... (7) 

r— »oo 

In IA.2I it is argued that taking rj = corresponds to a shift in the parameters of the 
displacement operators (with some constraints on the choice of new parameters), such that 
we may rearrange the \r])s regularisation to show that it approaches the NOPA regularisation, 
with tanhr = For = we therefore have: 

|r? = 0), = iV2e^'''''|00). (8) 

2.2 Tripartite states 

A suitable analogue of equations ([1]) or ([H]) which has the features required of an entangled 
state, which we analyse in detail below, is defined by: 

with normalisation \N^\'^ = \{s^ — 1)^ (s^ — 4) j-'^/^/s^. For the case i] = rj' = rj" = 0, the 
tripartite EPR-like state becomes: 

\r, = rj' = T]" = 0)s = iV3e^^"''''+"'''+''''^|000). (10) 



Note here that, while the set of states (jlOp belong to the well known pure, fully symmetric 
three-mode Gaussian states, the more general case of ([9]) where the parameters rj, rj and rj' 
are retained is not symmetric, since the parameters can all differ. For discussion of Gaussian 
states in relation to entanglement in CV systems, see |21] and references therein. 

Whereas the bipartite state \r\)s was a simultaneous eigenstate of {X\—X2) and (^1+^2)5 
in the tripartite case the choice of relative variables is no longer immediately apparent. In a 
similar manner to the derivation of ([1]), it is readily established using manipulations of the 
type: 

ae^ = |a - {A, a\ + \{A, [A, a]] + . . .| , (11) 
that generically \r],'q\rj")s is an eigenstate of the following combinations: 



(c-l(at + 6t))|r/,r/',/). = Vh, V, r?"),. (12) 

From this it is clear that different values of s will dictate limiting cases wherein \rj,rj' ,r]")s 
becomes a singular eigenvalue of various choices of relative variables. (Note that in the bi- 
partite case we could have introduced \r])s as \r],r]')s analogously, recovering ^ in the case 
r]' = —rj*.) Keeping s general, the eigenvalue equations become: 



Is: 



l^(^s + -^{X,-X2) + ^[s--^{P,-P2)M,v'')s = {V-V')hv',ri" 
±(^s + ^yX2-Xs) + ^(^s-^yP2-P3)hv',v")s = iv'-r]")hv',v")s, 
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From (I13|) it is clear that the singular cases will occur for s = 1 and s = For the case 
s = 1 we evidently have a singular eigenstate of the relative coordinates, while remaining a 
squeezed state of the total momentum. Conversely, for s = \/2 we have a singular eigenstate 
of the total momentum, but a squeezed state of the relative coordinates. 

If we construct mode operators corresponding to the Jacobi relative variables and the 
canonical centre-of-mass variables, say 

arel = ^(Xi-X3) + ^(Pl-P3), 

Wei = ^(Xi + X3-2X2) + ^(Pl+P3-2P2), 

acm = (Xi + X2 + X3) + -J= (Pi + P2 + Ps) , (14) 
then we find, from ()13p for general s: 



sarei + -al^i) \ri,ri ,7] }s = -j= [i] - i] ) \ri,ri ,ri }s, 



(^sbrel + ^^L) I^' v', v")s = (?7 - 2r/' + rj") |r/, rj' , 7]").,, 



_ 2 ^ 
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\r],v',v")s = -j={r] + r]' + r]")\rj,r]',r]")„ (15) 



from which it is again obvious that for s = 1 or s = \/2, canonical combinations arise in 
the first two, and last cases respectively. On the other hand, the non-canonical combinations 
appearing for s = 1 in the third, and s = \/2 in the first two cases, indicate that the squeezing 
parameters have the values ^ In 3 in each instance. 

Having established the structure of the tripartite EPR-like states, we can examine their 
behaviour when applied to Wigner functions, and the consequences of using these states in 
CHSH inequalities. 



3 Tripartite Wigner Function 

The Wigner function [221 [23] , was an attempt to provide the Schrodinger wavefunction with 
a probability in phase space. The time-independent function for one pair of x and p variables 
is: 

W{x, p) = — dy^* (x + y)4){x - y)e^'Py/^. (16) 
Trn J -00 

Alternatively, it has been shown [24[ [25] that a useful expression of the Wigner function is 
in the form of quantum expectation values. For modes, the Wigner function for a state {ip) 
may be expressed as the expectation value of the displaced parity operator, where the parity 
operator itself performs reflections about phase-space points (aj), where aj = -^{xj +ipj), 
with j = 1,2, . . . , N denoting the mode, and: 

W{ai,a2,...,aN) = (f )^ (n(ai, 02, • • • , OAr))- (17) 

The displaced parity operator is: 

n(ai,a2,...,ajv) = <E)f=iD,{aj)i-irw]{a,), (18) 

where rij are the number operators, and for each mode the Glauber displacement operators 
are of the form: 

D{a) = e°'^^-°*«. (19) 
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When we express the Wigner function in the form of (I17p . we can derive a set of these 
functions to construct the inequahties that wih be discussed in section [H The tripartite 
Wigner function for \r],r]' ,r]")s becomes: 

X (000| exp Q {ri*a + r]'*b + ri"*c) +^{ab + ac + bc)j 

^ ^aa^ -a* a ^fib^ - fS* b ^-yc^ --y* c ^_-^^na+ni,+nc ^a* a~aa^ ^I3*b-I3b^ ^-y* c--yc^ 

X exp Q (^r^at + 7?'6t + r]"c'<^ + ^ (^a+ftt + a^^c^ + ft+ct^^ |000). (20) 

We evaluate such matrix elements by commuting mode operators with the parity operator 
and rearranging using BCH identities before casting the operators into anti-normal ordered 
form. Then a complete set of coherent states is inserted and integrated over. As indicated 
in section [XT] and shown in \A.2\ we can absorb the rj, rj', ry" parameters by shifting the dis- 
placement parameters up to a factor: Wri^n\n"{a' , (3' ,'y') = -E(a', /?', 7', ry, ry', r/")Wo,o,o(«i 7)- 
We are free to choose instances where E{a' , (3' ,r],r]' ,rj") = 1, and henceforth we assume 
rj = Tj' = r]" = unless otherwise stated, and write simply VF(a,/3, 7). With this shift in 
displacements, the Wigner function for our tripartite state becomes: 

W{a,f3,j) = ^exp(^-^-J^-^[C7i(|ap + |/?P + |7|^) 

+ C2(a/3 + 07 + /37 + a* (3* + a*7* + /3*7*) 
+ C3{af3* + 07* + (3a* + f3^* + 7a* + 7/?*) 

+ C4(a2 + /J^ + 7^ + a*^ + + 7*^)] ) , (21) 



where 



Ci = -2(/-s^-4), C2 = 4s2(s4_2), C3 = -4s\ d = 



The most important point to note here is the emergence of mixed conjugate/non-conjugate 
pairs, which do not appear in the Wigner function for the second-quantised NOPA-like optical 
analogue (see lA.ip . To make the behaviour of the Wigner function in the asymptotic region 
clearer, the parameters a, /? and 7 are written in polar form, a = lale**^" etc. The Wigner 
function thus becomes: 



+ 2C2 {\a\\(3\ COs{(j)a + 0/3) + COs{<j)/3 + <j)^) + |7||a| COs(0-y + (j)a)) 

+ 2C3 {\a\\f3\ COs{<j)/3 - (pa) + COs((/)^ - (j)i3) + |7||a| cos(0^ - (pa)) 

+ 2C4 (lap cos(2</.„) + cos(20^) + \jf cos(20^))] ^ . (22) 

4 CHSH inequalities and violations 

The CHSH inequalities [6] are CV generalisations of the original Bell inequalities which were 
set up to test the scheme proposed by Einstein, Podolsky and Rosen (EPR) in 1935 [1]. Con- 
sidering the measurement of an entangled pair of particles performed after they have been 
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separated such that no classical communication channels are open when the wavefunction 
collapses, EPR posited that either quantum mechanics must be incomplete, with room for a 
hidden variable theory, or spatiotemporal locality is violated. Bell showed that hidden vari- 
ables were not permitted if we preserve both the assumptions of accepted theory - specifically 
locality - and the probabilities predicted by quantum mechanics. 

Generalised A^-mode Bell inequalities - CHSH inequalities, in terms of the Bell operator 
expectation values Bjsf - exist (see for example [16J). In their bi- and tripartite form we can 
apply these to our regularised EPR-like states. Following [inittS], the CHSH inequalities for 
the bi- and tripartite forms are the possible combinations: 

B2 = n(o,o) + n(o,/3) + n(Q,o) -n(a,/3), (23) 

1^21 < 2, 

B3 = n(0,0,7)+n(0,/3,0)+n(a,0,0)-n(a,/?,7), (24) 
l^sl < 2. 

In [13], the CHSH inequality constructed with the tripartite NOPA-like state (see equation 
()34p ) is maximised by taking an all-imaginary substitution a = /3 = 7 = i^/j, where J is some 
distance measure. In the bipartite case (Figure [1]), if we look in the region s — > 1"*" for ([8|), 
an all-imaginary substitution obviously gives exactly the same maximum violation as NOPA 
with r 00, both having a maximum value of i?2^"^ ~ 2.19 \10\ I13j. The figure shows clearly 
that the value of B2 increases as s — > 1"*" and J ^ 0. 

In the tripartite case however, we must examine the wealth of other possible choices which 
extremise the inequality. From the form of the Wigner function in (j22|) . however, there are 
some clear choices that will minimise the last term in B^. Choosing all the phases <pa — 
(pfB = (p-y = f ) ^-nd all the magnitudes \a\ = \j3\ = \j\ = such that all the parameters are 
imaginary {i^/j), equation ([22]) becomes: 

Ty(i^/J,0,0) = Ty(0,iVj,0)=iy(0,0,iVj) = ^exp(^-^^^^^y^^^, 

W{i^,i^,i^) = Aexp(- ^y_+^^ y (25) 

Consequently B^ is (from (pill ): 

r J{s^-s'^ + 2) 1 r 3J(s2 + 2)l 

In the region s ^ 1"^, -B3 never reaches a value greater than 2 (Figure [2]). A violation 
corresponding to the EPR limit s ^ l"*" can be found by making the choice a = —(3 = —\fj\ 
7 = 0, for which (j22p gives: 

H^(-\/j,0,0) = W^(0, 0) = exp 




vf(-Vj,Vj,o) = 

T^(0,0,0) = 1, (27) 
and Bj, becomes (Figure [3]): 

As s ^ 1"^, J ^ 0, the maximum value is B^""^ ~ 2.09, which can be checked both analytically 
and numerically. 
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However, what is more interesting still is exploring an auxiliary regime of the regulator, 
s — > \/2+, in equation ()26|) . This is shown in Figure [H Analytically, we can approximate the 
maximum to the lowest order in s — \/^ = e by writing = 3x — x^, where x = exp(— 4J/3e^), 
and A = 9, with maximum 

^max ^ (X-l) (^^^ ^ 2.32. (29) 

.„ 1 . — 
at X = yj) . This can be confirmed numerically for s v 2+, J ^ 0. The values of B^"'^ 

correspond exactly to those calculated for the experimentally verified NOPA-like states, whose 

maximisation as r ^ oo is also governed by (I29p . 



5 Discussion 



In this paper we have analysed tripartite CV entangled states which are natural generalisations 
of the classic bipartite EPR-type states (for two systems with canonical variables Xi, Pi, X2, 
P2). Given the necessity of working with normalisable states which still approximate the ideal 
EPR-type limit for practical implementation of CHSH inequalities, we examined a family 
of such regulated states parameterized by a regulating parameter s. This family of states 
was compared with those relating to multipartite NOPA-like states. The NOPA states have 
been shown to manifest CHSH violations, and have the advantage of being directly accessible 
by experiment via standard quantum optics protocols such as multiparametric heterodyne 
detection techniques and beam splitter operations. However, as an extension of a direct 
transcription of the EPR paradox, this new family of regularised states provides an alternative, 
systematic description of the approach to the ideal EPR states for relative variables. 

By finding expressions for the eigenstates of the regularised tripartite CV EPR-like states 
it became apparent that there are two regimes of the regularisation parameter in which these 
states become singular: in one case (s — > 1) we have a singular eigenstate of the relative 
coordinates while remaining squeezed in the total momentum; in the other, s — > \/2 limit we 
have a singular eigenstate of the total momentum, but squeezed in the relative coordinates. 
In these two regimes we have explored CHSH inequalities via Wigner functions regarded as 
expectation values of displaced parity operators. Violations of the tripartite CHSH bound 
(-B3 < 2) are established analytically and numerically, with B3 = 2.09 in the canonical regime 
(s — > 1"^), as well as B3 = 2.32 in the auxiliary regime (s \/2+). 

Related tripartite entangled states have recently been constructed by Fan [26]. Although 
these states are also accessible by standard quantum optics techniques, they are not true 
generalisations of 'EPR' states. In this case, while they diagonalise one centre-of-mass variable 
(for example, X1+X2 + X3), they are coherent states |27] of the remaining relative Jacobi 
observables (that is, they diagonalise their annihilation mode operators a, b in contrast to the 
s — > \/2 limit of our EPR-type tripartite states, which as stated above turn out to be squeezed 
states of these relative degrees of freedom (eigenstates of a linear combination -^(20 + a^), 

-^(2b -|- b^) in the relative mode operators, with the value ^ ln3 for the squeezing parameter). 
In the case of the tripartite entangled states of [2B], no regularisation has been given. A 
construction of true multipartite ideal EPR states has, however, been provided in [181 I19j . 
with a second-quantised form for the tripartite state (|48|) . which may be compared with the 
form for the NOPA-like state ()33p . Although the Wigner functions for the tripartite NOPA- 
like states show peaks at zeroes of Xi — Xj and Pi + Pj for all distinct pairs i,j which 
does not appear to be consistent with simultaneous diagonalisation of commuting observables, 
it can be inferred from the agreement of (j48p and (j33p that indeed in the infinite squeezing 
limit, tanh(r) = 1, and with relative parameters equal to zero, the NOPA-like state does again 
tend to the ideal EPR state. 
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Since the NOPA-like states are constructed with a view to experimental realisabihty, and, 
in the bipartite case, to manufacture the specific properties of the Wigner function, this new 
suggestion for a regularisation stemming from a direct transcription of the EPR paradox in 
terms of the simultaneous diagonalisation of commuting observables could be seen as a more 
general or fundamental description. It also considers in more detail the specific instance of 
tripartite EPR-type states, compared to the comprehensive [19] which finds n-partite rep- 
resentations of entangled states through their Gaussian-form completeness relation without 
exploring regularisations and Wigner function properties. As the proposed EPR-type regu- 
larised state produces a different Wigner function from the NOPA-type, with two singular 
limits, this paper's proposed regularisation may potentially suggest that alternative experi- 
mental ways to achieve the violations of the CHSH inequalities are possible. For a review of 
Gaussian states, and discussions of the realisabihty of entangled states, we refer to [2T | l28 l [29] . 
and references therein. It will be worth investigating the full extent of the constraints placed 
on the choices of displacement parameters entailed by the shift in r] (see IA.2p . The current 
discussion might also easily be extended to include a presentation of the alternative bipartite 
starting point of conjugate variable choice Xi + X2 and Pi — P2. [17j discusses the canonical 
combinations for any number of modes, but in our case it is reasonable to assume that an 

1 J J 



A^-partite generalisation would be of the form exp y-^ \ Yli<j 'H^jjj ■ would also expect 
that these states would admit standard completeness relations in the singular cases. 

In conclusion, we have presented a rigorous extension of Fan and Klauder's general EPR- 
like states to the regularised tripartite CV case for relative variables, and highlighted the 
connection to current quantum optics implementations. The CHSH inequalities constructed 
with component Wigner functions for this case show significant violation of the classical bound, 
and the different choices of regularisation parameter making the state singular illustrate an 
interesting new feature of the structure of generalised CV EPR-like states. 
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Deriving the tripartite NOPA and EPR-like Wigner functions 

In this appendix the tripartite NOPA and \rj,rj' ,r]")s Wigner functions are derived for com- 
parison. 

A.l Second-quantised form and Wigner function of tripartite NOPA-like 
state 

Applying two phase-free beamsphtters at specified angles acting on one momentum squeezed 
state and two position squeezed states of mode 1, 2 and 3 respectively, states that the 
tripartite NOPA-like states can be derived from the following expression: 



\NOPA^^'^) 



^23 (^) Bi2 ^arccos-^ 



X exp [ - I a" — a 



t2 



exp ( — I 



exp 



.t2 



1 000) .(30) 



We can therefore use the following formula quoted in |3lj0 for the squeezing operator S{z) 
(where z = e*^) with the Baker-Campbell-Hausdorff (BCH) relation: 



S{z) 



exp 



exp 



-iza 



e*^ tanh r)a 



t2 



t2 



exp 



-2(ln coshr)(-a^a-|--' 



exp 



1 



(e"*^ tanh r)a' 



(31) 



The following beamsplitter operatiorH can then be applied, where 9 here refers to the 
angles 7r/4 and arccos(l/\/3) for the i?23 and B12 splitters respectively: 

T3 ia\ f a ^ acos0 + bs'me , , 

[ — > — asmt' -1- ocos ^ 

and normalising, the tripartite NOPA is expressible in second-quantised form as: 
\NOPA^^'^) = (l-tanh2(r))^/^ 

X exp f-^tanhr {a^"^ +h'^'^ + 0^"^^ ^tanhr [h^ + a) + a) c^yj |000i(33) 



Using this state to derive the Wigner function of the form (jlTD . the complex exponential 
that is produced may be rearranged using the formula (03]). The tripartite NOPA Wigner 
function then becomes: 



Wnopa 



+ 



2Y 

- exp 
nj 

2 tanh r 

3 (1 — tanh ^r) 

8 tanh r 

3(1 — tanh ^r] 



1 — tanh 2? 



(a*2 + ^*2 ^ + ^2 ^ ^2 ^ ^2^ 



(a/3 + /37 + 7a + a* (3* + ^7* + 7*"*) 



^exp{(-2cosh(2r))(|a|2 + |/3|2 + |7|2) 

sinh(2r) (a*2 + p*^ + 7*2 + „2 ^ ^2 ^ ^2^ 

4 

-f-- sinh(2r) (a/5 + p-f + ja + a* (3* + (3*-f* + 7*0*) 



(34) 



■^Note the misprint in the sign of the last exponential in [31]; see |32) 

^An additional overall relative sign (180° phase shift) between the two modes has been omitted; see for 



example [33] 
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This is the result quoted in [13], and further exphcation can be found in that paper. This 
function should be compared with the Wigner function of our regularised tripartite EPR-like 
state, \r],r]' ,r]")s. Further details of that derivation are given below. 



A. 2 Derivation of conditions for r] = 

In the interest of brevity, the conditions for shifting the rj parameters are shown below for the 
biparite case. However, the analysis extends in an obvious way to the triparite case. After 
BCH and anti-normal ordering, the bipartite Wigner function (tripartite given in equation 
(f^n|) ) becomes: 

W{a,f3) = ,(7?,7? le'l"l'e^l^l'e-^°*"e-'^*''e'°"^'^''^X -l)"-+"''|r?,V), 

F(a,/3) 

= (00| exp (-^|r/P - + + ^r]'*b + ^ab^ F{a, (3) 

X exp - ^\vf - -m^ - -/b^ + ^a^b^) |00). (35) 

We then make the generic substitutions 

a = a' + A{r],s), 

P = f3' + B{rj',s), (36) 

into F{a, 13). To find the expressions for A{r], s) and B{r]' , s) that will allow us to set r] = rj'{= 
r]") = 0, we solve the following: 

2a'A* + 2a*A + 2\A\^-2A*a + 2Aa) = ■^\r]\^ - -ji* a + -j]a) , 



s s 



2P'B* + 2P'*B + 2\B\'^ -2B*b + 2Bb^ = -XrW]^ - -r]'*b + -r]'bl (37) 

s s 

Thus we can see that, if we allow the constraints 

a ri an 

— - + = 0, 

iil+llk = 0, (38) 

s s 

(i.e. a' real and rj imaginary or vice versa), then the expressions for A(r],s) and B(r]',s) 
become 

A{v,s) = ^^, A*{rj,s) = ^^, 

B{r,',s) = ^^, B*{r^',s) = ^. (39) 

Therefore, up to a factor, taking r] = r]'{= r/") = corresponds to a shift in the parameters 
of the displacement operators a' = a — ^, /?' = /? — ^ (and 7' = 7 — ^ in the tripartite case). 
For the tripartite Wigner function as used in section ([3|), this can be expressed as: 

Wrj,r,',v"i»'^f3'^i) = E{a',l3',^',rj,rj',r,")Wo,o,o{»,(3,j) (40) 
E{a' , P' ,ri,ri' ,ri") = exp (a'r/* + a'* 7/ + /3'r/'* + + 7'??"* + 7'* 77") j (41) 
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A. 3 Details of derivation of tripartite \ri,ri' ,'r]")s Wigner function 

The second-quantised EPR-like state is expressed as ®. This is used to find the Wigner 
function in the form of ()20|) . By commuting mode operators with the parity operator and 
rearranging using BCH identities, the expression becomes, in anti-normal ordered form: 

3 



W 



X (000| exp Q {r]*a + r]'*b + r/"*c) + ^{ab + ac + bc)^ 

^ g2|a|2g2|/3|2g2|7|2g-2a*ag-2/3*feg-27*Cg2aatg2/36tg27ct 

X exp (^ria^ + rj'b^ + 7?"c"f) + ^ (^a^b^ + a^c' + b^c^^^ |000). (42) 



In anti-normal ordered form we may insert a complete set of coherent states J \ u,v, w) (n, v, 
such that we may rearrange the exponential according to the formula [3; 



w 



d?ud?vd?w 



n 



IT 



exp 



A B 
C D 



z 
z 



det 



det 



C D 

A B 

C D 

A B 



exp 



2^^'^''*\ C D 



A B 



exp 



C D 
A B 



z 
z* 

V* 



(43) 



where matrices A and D must be symmetrical, and C = B^ . In this instance 

(z,z*) = {u,v,w,u* ,V* jW*) , 



and we have 



1 



1 



-7] 



C D 
A B 



'* 


-2(5* 


\"* 


-27* 


1 


ri + 2a, 


1 , 

— V 4 






s 




s 




s 


/ 


1 











1 









1 





1 

"? 





1 










1 




1 

~^ 










1 


1 

"7 


1 










1 

"7 










1 





\ 




1 











1 / 



(45) 



with inverse 



C D 
A B 



-1 



Note also that 



det 



C D 
A B 



,12 





1 


1 


2s-2 


^4-2 


s4-2 \ 


1 




1 


s4-2 


2s-2 




1 


1 




34-2 


s4-2 


2s-2 


2s-2 


s4-2 


^4-2 


5^-3 


1 


1 


s4-2 


2s-2 


s4-2 


1 


s4-3 


1 




s4-2 


2 5-2 


1 


1 


s^-3 j 



es^ + Os"^ -4)/s 



121 



(46) 



(47) 
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such that the cancel in the Wigner function. 

Prom the argument in IA.21 we assume that rj = rj' = rj" = unless otherwise specified, 
and continue to use W{a, f3,'j). This gives equation (|21|) . which may now easily be compared 
with the Wigner function derived for the NOPA-like case (equation (I34p ). Further discussion 
of similar manipulations of the Wigner function can be found in [36] , 



A. 4 Tripartite entangled state from 

Equation (27) in [19j provides the ideal EPR state for the tripartite entangled state: 



4 



^ + ^ E «1 + ^(4 - 24 + 4) + 4^(4 + 4 - 24) + s^] m, 



3 .=1 



A ^ -y-^(ei + el-6e3), 



5 



2 3 .3 



E (48) 



6 

i<j=l i=l 
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Figures 




Figure 1: Plot of bipartite s-modified CHSH, with an all-imaginary choice for a and /?. Reaches 
a maximum value of 2.19 as s — >^ 1 and J ^ 0. Note that this is equivalent to the NOPA 
case. 




Figure 2: Tripartite s-modified CHSH. With an all-imaginary choice for a, P and 7, -B3 never 
reaches a value greater than 2 as s ^ 1. 
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